function test()
A=[1,4,5;4,18,26;3,16,30];
[L,U,P]=LU_factorization(A);
display('input:');
display(A);
display('output:');
display(L);
display(U);
display(P);
end